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We have carried out systematic Macroscopic Quantum Tunneling (MQT) experiments on Nb/Al- 
A10cc/Nb Josephson junctions (JJs) of different areas. Employing on-chip lumped element inductors, 
we have decoupled the JJs from their environmental line impedances at the frequencies relevant for 
MQT. This allowed us to study the crossover from the thermal to the quantum regime in the low 
damping limit. A clear reduction of the crossover temperature with increasing J J size is observed 
and found to be in excellent agreement with theory. All junctions were realized on the same chip 
and were thoroughly characterized before the quantum measurements. 

PACS numbers: 74.50.+r, 85.25.Cp, 74.78.Na 



I. INTRODUCTION 

Since the gauge- invariant phase over a Josephson junc- 
tion (JJ) if is a macroscopic variable, circuits containing 
JJs have been used as model systems for the investiga- 
tion of quantum dynamics on a macroscopic scale. This 
research has recently led to the development of differ- 
ent types of superconducting quantum bits^"— , which are 
promising candidates for the implementation of quantum 
computers. The starting point of this field was the ob- 
servation of Macroscopic Quantum Tunneling (MQT) in 
Josephson junctions in the 1980s^. In such experiments, 
the macroscopic variable ip is trapped in the local min- 
imum of a tilted washboard potential, before it tunnels 
through the potential barrier and starts rolling down the 
sloped potential. Since this running state is equivalent 
to the occurrence of a voltage drop over the junction, 
such tunneling events can be experimentally detected. 
MQT — often referred to as secondary quantum effect — 
is the manifestation of the quantum mechanical behavior 
of a single macroscopic degree of freedom in a complex 
quantum system. Furthermore, it is the main effect on 
which all quantum devices operated in the phase regime 
(such as phase qubits and flux qubits) are based. Conse- 
quently, the detailed understanding of MQT is not only 
interesting by itself, but also important for current re- 
search on superconducting qubits operated in the phase 
regime. In this article, we report on a systematic ex- 
perimental study of the dependence of the macroscopic 
quantum tunneling rate on the Josephson junction area, 
which to our knowledge has never been performed be- 
fore. As usual 6 8 , we measure the rate at which the es- 
cape of ip out of the local minimum of the washboard 
potential occurs as a function of temperature. At high 
temperatures, the escape is driven by thermal fluctua- 
tion over the barrier while it is dominated by tunneling 
at low temperatures. This leads to a characteristic sat- 
uration of the temperature dependent tunneling rate be- 



low a crossover temperature T cr , which is the hallmark of 
MQT. The rates above and below crossover are affected 
by the dissipative coupling to the environment of the J J, 
which is commonly accounted for by a quality factor Q in 
theoretical descriptions. A major goal of the presented 
study was to keep the influence of Q on the rate constant 
while varying the junction area, so that a change in the 
observed escape rates could be clearly assigned to the 
changed J J size. For this purpose, we work in the under- 
damped regime of large Q, which is only possible if the J J 
is to some extent decoupled from its low-impedance en- 
vironment (i.e. the transmission line leading to the JJ). 
We achieve this by employing on-chip lumped element 
inductors. 

This article is organized as follows: First, the physical 
model of MQT is discussed and the theoretical expec- 
tations for varying junction size are given. Second, the 
procedure and setup of measurement are described. Af- 
terwards, the investigated Josephson junctions are char- 
acterized carefully, and finally, the results of the MQT 
measurements are presented and discussed. 



II. MODEL AND MACROSCOPIC QUANTUM 
TUNNELING 



A. General Model 

The dynamics of a JJ is usually described by the RCSJ 
(resistively and capacitively shunted junction) model^iiS. 
The current flowing into the connecting leads comprises 
in addition to the Josephson current Ij = I c sin ip (I c de- 
notes the critical current of the junction) a displacement 
current due to a shunting capacitance C and a dissipative 
component due to a frequency dependent shunting resis- 
tance R. For a complete description, the electromagnetic 
environment given by the measurement setup can be in- 
cluded in the model parameters. In our case, R will be 
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influenced by the environmental impedance while C can 
be regarded as solely determined by the plate capacitor 
geometry of the J J itself. In any case, the bias current / 
is composed of 



I = I c sin (f 



1*2 • 

R2^ 



(i) 



where cp is the gauge-invariant phase difference across 
the junction and <3>o = h/2e is the magnetic flux quan- 
tum. The dynamics of cp as expressed by ([1]) is equally 
described by the well-studied Langevin equation 



M(p + r]M(p + ^ = £(t) . 
dip 



(2) 



which describes a particle of mass M = C(^o/2tt) 2 in a 
tilted washboard potential 



U ((p) = Ej (1 - cos p - j(p) 



(3) 



exposed to damping 77 = 1 / i?C and under the influence of 
a fluctuating force £(£). The strength of is linked to 
temperature and damping by the fluctuation-dissipation 
theorem. Furthermore, 7 = I/I c denotes the normalized 
bias current while Ej = <&qI c /2k is called the Josephson 
coupling energy. For 7 < 1, if thermal and quantum 
fluctuations are ignored, the particle is trapped behind a 
potential barrier 



AU = 2Ej v 1 - 7 



7 arccos 7 



(4) 



and the JJ stays in the zero-voltage state. In the po- 
tential well, the phase oscillates with the bias current 
dependent plasma frequency 



Up u p o (l - 7 2 ) 



2\V4 



2\l/4 



(5) 



To complete the list of important system parameters 
given in this section, we introduce the quality factor 
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which is conventionally used to quantify the damping in 
the JJ. 

At finite temperatures, the thermal energy ksT (ks 
being Boltzmann's constant) described by in ((2]) can 
lift the phase particle over the potential barrier before the 
critical current 7 = 1 is reached, so that the particle will 
start rolling down the potential. This is called prema- 
ture switching and the observed maximal supercurrent 
7 SW < I c is called the switching current. When the phase 
particle is rolling, the JJ is in the voltage state, since 
a voltage drop according to (p = (2tt/^o)V is observed. 
The thermal escape from the potential well occurs with 
a rate— ^ 
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where a t is a temperature and damping dependent pref- 
actor, which will be discussed in more detail in Sec. Ill Ci 
For T —> 0, where T t h —> 0, premature switching will 
still be present due to quantum tunneling through the 
potential barrier. As the phase difference over the J J is a 
macroscopic variable, this phenomenon is often referred 
to as "Macroscopic Quantum Tunneling" (MQT). This 
means that by measuring the switching events of a JJ 
for decreasing temperature, one will see a temperature 
dependent behavior (dominated by the Arrhenius factor 
exp(— ACI/fo^T) in ©) until a crossover to the quantum 
regime is observed. The crossover temperature T cr is ap- 
proximately given byii^ 
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where h is Planck's constant. We can write the quan- 
tum tunneling rate for temperatures well below crossover 



as 
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(9) 



where a q = y/%§4nMJ/hw p exp(1.430/Q) and B = 
(36AU/5hujp)(l + 0.87/Q). In the limit of large Q, the 
escape rate is expected to approach the temperature in- 
dependent expression ([9]) quickly 15 -^ once the tempera- 
ture falls below T cr . The rates in ([7]) and (j9j) are func- 
tions of the normalized bias current 7 via @ and (|5j). 
The crossover to the quantum regime can be nicely visu- 
alized by measuring the bias current dependence of the 
escape rate T(I) for a sequence of falling temperatures. 
The data are then described over the whole temperature 
range by the thermal rate ([7]) with the temperature as a 
fitting parameter. In this way, one obtains a virtual "es- 
cape temperature" T esc , which can be compared to the 
actual bath temperature T. In the thermal regime, one 
should obtain T esc = T while in the quantum regime, one 
should get T esc = T cr = const. 



B. Influence of JJ Size on MQT 

The crucial element of a Nb-based Josephson junc- 
tion as employed in this work is the Nb/Al-AlO^/Nb 
trilayer. For a JJ, I c = j c • A and C = c • A : where the 
critical current density j c and the specific capacitance 
c are constant for a given trilayer and A is the area of 
the junction. Hence, by reformulating (|5j), we find that 
Upo = Y // 27rj c /<l>oc, meaning that for JJs fabricated with 
the same trilayer, the plasma frequency does not depend 
on their size. So at first sight, the crossover temperature 
JHJ) should also be independent of the JJ size. In real- 
ity, however, the problem is more subtle, as one needs 
to take into account at which normalized bias current 
7 cr the quantum tunneling rate ((9]) becomes significant. 
Since the height of the potential barrier AU oc Ej oc A 
is proportional to the JJ size, a significant tunneling rate 
should be reached at different 7 cr values for junctions of 
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0.95 0.96 0.97 0.98 0.99 
7=111 



1.00 



TABLE I. Parameters of the investigated samples and cal- 
culated 7cr and T cr values. The critical current density ac- 
counts for j c w 650 A/cm 2 while the specific capacitance of 
c — 55 fF/um 2 was recently determined by the measurement 
of Fiske steps on a trilayer prepared in our lab under identical 
conditions. Parasitic capacitances due to idle regions next to 
the JJs have been taken into account for all calculations in 
this article. For the calculation of T cr , the actual measured 
critical currents were used. 



Sample 


Diameter d (um) 


7cr 


T cr (mK) 


Bl 


1.9 


0.965 


371 


B2 


2.55 


0.977 


323 


B3 


3.6 


0.988 


291 


B4 


3.8 


0.988 


277 



FIG. 1. Theoretically calculated switching current distribu- 
tions P(I) (solid curves) and quantum tunneling rates T q 
(dashed curves) for samples Bl (left) to B3 (right) having 
different diameters d (for parameters see Tab. UJ) . The dif- 
ference in 7 cr (maximum position of P(I)), where quantum 
tunneling leads to escape from the potential well, is signifi- 
cant. The tunneling rates at these points (marked by arrows) 
are of the order of ~ 100 kHz for all samples. 



different size. These points can be estimated by theo- 
retically calculating (j9j and converting it into a switch- 
ing current histogram, as it would be observed in a real 
experiment. The probability distributions of switching 
currents P("f) can be obtained from the quantum rate by 
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where d'y/dt — (1/I C ) • (dl/dt) is a constant for the lin- 
ear current ramp chosen in our experiment. For the pa- 
rameters of the junctions investigated in this work (see 
Tab. UJ, the switching current distributions P(j) were 
determined with a quality factor of Q = 100 and a cur- 
rent ramp rate of 100 Hz, according to our experiments 
(see below). They are shown in Fig. [TJ where it can be 
seen that the 7 cr values (the positions of the maxima 
of the distributions) significantly and systematically in- 
crease with the junction size. Evaluation of the r g (7 cr ) 
values for the samples indicates that quantum tunnel- 
ing will be experimentally observable at a rate of around 
T q w 10 5 Hz. 

Subsequently, the expected crossover temperature was 
calculated from (|5J with 7 = 7 cr . The sample parame- 
ters as well as the expected 7 cr and T cr values are given 
in Table Ul It can be seen that due to the term in paren- 
thesis on the right hand side of (|8j, the crossover tem- 
perature systematically decreases for increasing junction 
size. The change in T cr is large enough to be observed 
experimentally. However, such a systematic study of the 
size-dependence of T cr has never been carried out before. 



C. Influence of Damping on MQT 



The quality parameter Q is frequently employed to 
describe the strength of the hysteresis in the current- 
voltage characteristics of a JJ. In this case, one often 
takes Q = uj p oR sg C with R sg being the subgap resis- 
tance of the junction. Here, Q is size-independent, as 
R sg oc 1/A and C oc A. In the context of MQT, however, 
the dynamics takes place at a frequency of uo pi so that 
a complex impedance at that frequency Z(uo p ) has to be 
considered. For an MQT experiment, where the phase 
and not the charge is the well-defined quantum variable, 
the admittance Y(uj p ) will be responsible for damping^, 
so that R in (|6j will be given by R = l/Re(F). 

If the junction was an isolated system, the value of R in 
the context of MQT would be determined by the intrin- 
sic damping in the zero- volt age state. The value which is 
typically taken as a measure for this is the maximal sub- 
gap resistance i? S g,max, which is simply the maximal re- 
sistance value which can be extracted from the nonlinear 
subgap branch of the current- voltage characteristics^^. 
In most experiments however, the electromagnetic envi- 
ronment of the JJ can be assumed to have an impedance 
that is real and accounts for Zq w 100 fi, correspond- 
ing to typical transmission lines 7 . As furthermore Zo <C 

R. 



sg,max 



and both contributions are in parallel (see Fig. 
[2^), we can simply write Q = uj p ZqC in this case. 

Evidently, for junctions having a small capacitance (as 
in our experiment), the quality factor Q = uo p ZqC will 
be limited to Q < 10 and additionally depend on the JJ 
size like C oc A. As we want to investigate the pure influ- 
ence of the JJ size on MQT, we would like to obtain very 
low damping as well as similar damping for all investi- 
gated junctions. In the implementation of phase qubits, 
current biased Josephson junctions have been inductively 
decoupled from their environment by the use of circuits 
containing lumped element inductors and an additional 
filter junction 3 . In order to keep our circuits simple, we 
attempted to reach a similar decoupling by only using 
on-chip lumped element inductors right in front of the 
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Metal powder filter 



FIG. 2. a) Typical impedance environment for switching 
experiments in a J J. As Zq <C R sg , the junction sees the 
impedance Zo at the plasma frequency, b) Lumped element 
inductors L can be used to decouple the JJ from the line 
impedance Zo, as discussed in the text, c) SEM micrograph 
of the electrode design used for the investigated junctions. 



JJs (see Fig. EJd). This setup leads to an admittance 



Y = 1/R S: 



l/(Z + iu>L). 



(11) 



As for (HQ), we find Re(Y) -> 1/R, 



sg,max 



in the limit 

uL oo, big enough lumped element inductances should 
decouple the JJ from the Zo environment and result in 
a high intrinsic quality factor Q = oj p R sg ^ max C even for 
switching experiments. Although it might be difficult to 
reach this limit in a real experiment, decoupling induc- 
tors should definitively help to increase the quality factor 
and move towards a JJ-size independent damping. 

The damping in the JJ influences the thermal escape 
rate (|ZJ) via the prefactor at < 1, which has been cal- 
culated for the first time by Kramers in 1940 1 — . In the 
limiting case Q —> (moderate to high damping), he 
found: 



a>t = «KMD 




1 

2Q J ' 2Q ' 



while in the opposite limit Q — >• oo (very low damping 
limit), he found: 



a>t — <^KLD 



36A/7 
5Qk B T 



More recently, Biittiker, Harris and Landauer— extended 
the very low damping limit to the regime of low to mod- 
erate damping finding the expression 23 



a t 



(Vl + 4/axLD + l) 2 



(12) 



Additionally, damping reduces the crossover tempera- 
ture according to 12 i 24 



• «KMD • 



(13) 
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Cu RF screen 
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Superconducting coil 



RF shield -± 

FIG. 3. Schematic overview of the measurement system. The 
superconducting coil and the sample are inside a magnetic 
shield consisting of three nested cylindrical beakers, the mid- 
dle one made from Pb, the two remaining ones from Cryop- 
erm. Furthermore, the entire dilution refrigerator is placed 
inside a /i-metal shield at room temperature. The 7r-symbols 
denote commercial 7r-filters. 



A possible way to determine the quality factor Q for 
such quantum measurements is to extract it from spec- 
troscopy data 7 . Unfortunately, for samples with such a 
high critical current density as used in our experiments 
described here, this turns out to be experimentally very 
hard. Hence, we will limit the analysis of the damping 
in our experiments to the MQT measurements. However, 
other groups have found a good agreement between the Q 
values determined by spectroscopy and by MQT— ^ and 
we hope to observe such a major increase in Q due to 
the decoupling inductors that minor uncertainties in Q 
should not play a role. 



III. SETUP AND PROCEDURE OF 
MEASUREMENT 

All samples were fabricated by a combined pho- 
tolithography / electron beam lithography process based 
on Nb/Al-AlO^/Nb trilayers. The trilayer deposition 
was optimized carefully in order to obtain stress-free Nb 
films. For the definition of the Josephson junctions, an 
Al hard mask is created employing electron beam lithog- 
raphy. This hard mask acts as an ideal etch stopper 
during the JJ patterning with reactive-ion-etching. Fur- 
thermore, it allows the usage of anodic oxidation even 
for small junctions, which would not be possible if a re- 
sist mask was used. After the anodic oxidation, the Al 
hard mask is removed by a wet etching process. Details 
of this Al hard mask technique and the entire fabrication 
process are discussed elsewhere 2 ^. 

Our measurement setup can be seen in Fig. [3l Spe- 
cial care has been taken in design of the filtering stages 
in order to reach a low- noise measurement environment. 
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800 mK 


o j = 
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FIG. 4. Top: The measured switching current histograms 
for sample B3 for selected temperatures. For increasing T, 
the switching currents decrease and the histograms broaden. 
Bottom: The plot obtained by applying (|15|) for the same 
sample. The fits allow to extract T esc as well as I c . 
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FIG. 5. Characterization of sample B3. The IV curve shows 
the high quality regarding the I c Rn ratio as well as low sub- 
gap currents. The inset shows a magnification of the subgap 
branch achieved by a voltage bias. The blue line illustrates 



how the value for R s < 



was determined. The fact that the 



current rises for decreasing voltage at V ~ 0.5 V is due to the 
fact that the junction jumps back to a supercurrent I ^ for 
V = 0. 



The goal of the measurement is to determine the escape 
rate T. In order to do so, we have measured the proba- 
bility distribution P(I) of switching currents. This was 
done by ramping up the bias current with a constant rate 
I = dl/dt and measuring the time t sw between 1 = 
and the switching to the voltage state with a Stanford 
Research 620 Counter, so that J sw = / • t sw could be 
calculated. An Agilent 33250A waveform generator was 
used to create a sawtooth voltage signal with a frequency 
of 100 Hz, which was converted into the bias current by 
a resistor of 47 kfi. In this way, for each temperature, 
I sw could be measured repeatedly. After doing so 20,000 
times, the switching current histograms P(I) with a cer- 
tain channel width AI were attained as shown in the 
upper part of Fig. 2J These histograms were then used 
to reconstruct the escape rate out of the potential well 
as a function of the bias current by employingi^ 

^-s'S <l4) 

With T at hand, we could now determine the escape 
temperature T esc by employing ([7]). In order to be able 
to rearrange this formula, we approximate the potential 
barrier in the limit 7 1 as AU = 4y/2/3- Ej- (1 -7) 3/2 , 
so that we find 

/ 2*T(I) \ 2/3 _ ( 4V2Ej \ V3 I c -I 

V n a t {I)uj p {I) ) \3k B T esc ) I c ■ [ > 

Hence, by plotting the left side of JTSJ) over the bias cur- 
rent /, we should obtain straight lines (see bottom part 
of Fig. 2]). Consequently, we can extract the theoretical 



critical current I c in the absence of any fluctuations as 
well as the escape temperature T esc by applying a linear 
fit with slope a and offset b. We then find 

T b 4>/2$o 

I c = -- and T esc = -=. . 

a QirkBCLvb 

Since I c enters ([T5]) via Ej and uo pl this fitting proce- 
dure has to be iteratively repeated until the value of I c 
converges. So strictly speaking, this procedure involves 
two fitting parameters, namely T esc and I c . However, 
it turns out that I c is temperature independent within 
the expected experimental uncertainty (for all our mea- 
surements, the fit values of I c vary over the entire tem- 
perature range with a standard deviation of only around 
0.09 %). Furthermore, the found I c values agree very well 
with the expected ones from the critical current density 
j c of the trilayer and the junction geometry. Altogether, 
it can be said that the results for the main fitting param- 
eter T esc should be very reliable. 



IV. SAMPLE CHARACTERIZATION 

The J Js were circular in shape and their geometries are 
given in Table H In order to characterize the samples, 
TV curves with current bias as well as IV curves with 
voltage bias were recorded (an example can be seen in 
Fig.[5|). The quality parameters for all samples are given 
in Table HI1 and indicate a very high quality. In the voltage 
bias measurements, two major current drops at voltages 
2 A/2 and 2 A/3 could be seen and attributed to Andreev 
reflections 26 . Below 2 A/3, we were able to extract values 
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of the maximal subgap resistance i? sg;max 
by the blue line in Fig. 



as illustrated 



TABLE II. Experimentally determined parameters for all in- 
vestigated JJs. The theoretical critical currents I c were ex- 
tracted from the MQT measurements. The i? sg ,max values 
were obtained as shown in the inset of Fig. [5ji. 



Sample 


Ic (liA) 


V gap (mV) 


I C R N (mV) 


-Rsg,max (k£~2) 


B-l 


19.1 


2.88 


1.75 


54.0 


B-2 


31.9 


2.88 


1.86 


73.0 


B-3 


68.1 


2.92 


1.93 


21.1 


B-4 


70.8 


2.90 


1.91 


31.5 



0.01 



0.008 ■ 



=3 0.006- 



0.004 ■ 



0.002 ■ 



Sample B1 

Sample B2 

Sample B3 

— Sample B4 



40 60 80 100 120 140 160 180 200 

Quality Factor Q 



V. RESULTS AND DISCUSSION 
A. Damping in the Junctions 

In order to decouple the JJs from their environmental 
impedance, the electrodes leading to the junctions were 
realized as lumped element inductors, as can be seen in 
Fig. [2t. This design was based on the layout that we 
recently used to successfully realize lumped element in- 
ductors for LC circuits in the GHz frequency range 27 . 
Furthermore, simulations with Sonnet 28 confirmed that 
the meandered electrodes indeed act as lumped element 
inductors at the relevant frequencies uj p (j cr ). The com- 
plex simulation with Sonnet gives an inductance of L/2 « 
1.65 nH (for one electrode) while the much simpler anal- 
ysis with FastHenry 29 yields L/2 w 1.8 nH. 

For each sample, the data were analyzed using a num- 
ber of different Q values in order to see if we could de- 
termine the experimentally observed damping. This was 
done by calculating the deviation of T esc from the bath 
temperature T in the thermal regime: 

AT 2 = Y, ( T esc-T) 2 (16) 

T>500 mK 

and finding its minimum value regarding Q. The corre- 
sponding values were then used for the sample analysis. 
It can be seen in Fig. [6] that the points of experimen- 
tally observed damping could be clearly identified. The 
evaluated Q values are given in Table IIII1 

In a preliminary experiment, we investigated MQT in a 
junction with a diameter of d = 1.9 um, a critical current 
of I c « 12 uA and low- inductance electrodes, which were 
simple wide lines and can be imagined as the envelope 
of the electrodes in Fig. [2t. We carried out a similar 
analysis to determine the damping and obtained a quality 
factor of Q = 4. Subsequently, we evaluated (j6j) and 
calculated an impedance of R = 99.8 ^, which is very 
close to the expected value of Z ~ 100 Q for typical 
transmission lines 7 . This means that with this simple 
preliminary design, the junction was in no way decoupled 
from the electromagnetic environment. 



FIG. 6. Determination of the experimentally observed quality 
factor Q for all samples. The curves are minimal when the de- 
termined Tesc values deviate the least from the corresponding 
bath temperatures T in the thermal regime T > 500 mK. 



The Q values obtained by using inductive electrodes 
(see Tab. HIT]) , however, show that we have drastically 
increased the quality factors with respect to the prelim- 
inary measurement. If we calculate the R values us- 
ing (J6j), we find that they are clearly above the typi- 
cal line impedance of Zq « 100 Q as well as the vac- 
uum impedance of 377 ^, which shows that we were in- 
deed able to inductively decouple the JJ from its usual 
impedance environment. As expected, the determined 
R values are still clearly below the subgap resistance 
max , indicating that we have not reached the limit 
ujL — >• oo. Instead, we are in the intermediate regime 
Zo <C R <C i? S g,max, leading to the fact that Q still ex- 
hibits a slight dependence on the J J size (see Table HI!]) . 
However, all JJs are in the low-damping regime, so that 
no influence of damping on the results should be present 
and differences in the experimental results should indeed 
be due to the JJ size. This can be seen by the fact that 
the damping related correction in T cr according to equa- 
tion ([T3]) is smaller than 1 % for all experimentally ob- 
served Q values. Altogether, we can state that we will 
be able to carry out our investigation of the size depen- 
dence of MQT with very low and nearly size-independent 
damping. 

In addition to the rather qualitative considerations 
above, we performed a quantitative analysis employing 



TABLE III. The experimentally determined values charac- 
terizing the damping for all samples. The L ca \ c values were 
determined from Q, the values given in Table Ull and equation 
(|11|) . They are in good agreement with the design value of 
L^3.3nH. 



Sample 


Q 


R (O) 


£ ca ic (nH) 


B-l 


76 


1571 


1.33 


B-2 


98 


1327 


1.39 


B-3 


132 


1015 


1.36 


B-4 


143 


1041 


1.43 
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Bath temperature T [K] 



FIG. 7. Calculated escape temperatures for sample B3 with 
and without an applied magnetic field. The inset shows the 
ic(<i>) modulation of this junction; the arrows indicate where 
the MQT data was obtained. For I m &gn = 31.6 mA, a clear re- 
duction of the observed crossover temperature T cr is observed. 



equation ([IT]) . If we use uo = co> p (7 cr ), take i? sg , m ax from 
Table [III and assume that Zq = 100 £1, we can calculate 
the decoupling inductance L CdL \ c for all samples. The val- 
ues, given in Table IHIi are a factor of around 2.3 — 2.5 
smaller than the simulation value of L w 3.3 nH, but of 
the right order of magnitude. For such a complex system, 
this is a surprisingly good agreement between simulation 
and theory on the one side and experimentally deter- 
mined values on the other side. In summary, we conclude 
that we have successfully demonstrated that decoupling 
of the Josephson junction from its environment is also 
possible using only lumped element inductors. 



B. Crossover to the Quantum Regime 

We now turn to the investigation of the crossover point 
from the thermal to the quantum regime and the influ- 
ence of J J size on it. As can be seen in Table HI we expect 
a clear reduction of T cr with increasing J J size. However, 
an experimental observation of lower crossover tempera- 
tures for smaller J Js having smaller critical currents could 
simply be due to current noise in our measurement setup. 
In order to exclude this, we artificially reduced the crit- 
ical current of sample B3 by applying a magnetic field 
in parallel to the junction area. While unwanted noise 
should now lead to an increase in the observed T cr , the 
physical expectation is a significantly reduced T cr due to 
the lower plasma frequency according to J5J). The result 
of this measurement can be seen in Fig. [71 and Table IIVI 
We found an agreement between calculated and observed 
crossover temperature down to T cr « 140 mK, which was 
the lowest temperature we examined. Hence, it is clear 
that we have a measurement setup exhibiting low noise, 
where the electronic temperature is indeed equal to the 
bath temperature. The lowest investigated temperature 
of 140 mK is clearly below any temperature needed for 




Bath temperature T [K] 

FIG. 8. Calculated escape temperatures for all samples. The 
crossover to the quantum regime is very clear in each mea- 
surement. The inset shows a magnification of the quantum 
regime. The reduction of the crossover temperature T cr with 
increasing J J size is clearly visible. 



the comparison of the JJs of different sizes with each 
other. 

Finally, we measured the switching histograms for the 
four JJs of different sizes and evaluated the escape tem- 
perature T esc and the theoretical critical current I c . This 
allowed us to determine the crossover temperature T cr 
and the normalized crossover current 7 cr = i sw ,cr/^c- We 
indeed found a clear dependence of the crossover tem- 
perature on the JJ size as can be seen in Fig. [8j To 
compare the experimental 7 cr and T cr values with the 
ones expected by theory, we now performed the theoreti- 
cal calculation described above using the experimentally 
determined Q values and equation ([T3]) . All experimen- 
tally determined values are in excellent agreement with 
theory, as can be seen in Table IIVI 

TABLE IV. The experimentally determined values character- 
izing the crossover from the thermal to the quantum regime 
in comparison with the theoretical expectations for all mea- 
surements. 



Sample 


W/c(o) 


^crjtheo 


7cr,exp 


^~cr,Q,theo 

(mK) 


7cr,exp 

(mK) 


B-l 


i 


0.965 


0.970 


368 


362 


B-2 


i 


0.977 


0.981 


322 


321 


B-3 


i 


0.988 


0.990 


290 


294 


B-3 


0.52 


0.984 


0.987 


223 


236 


B-3 


0.27 


0.979 


0.983 


172 


176 


B-3 


0.13 


0.973 


0.975 


129 


147 


B-4 


1 


0.988 


0.990 


276 


278 



VI. CONCLUSIONS 

We have carried out systematic Macroscopic Quantum 
Tunneling (MQT) experiments with varying Josephson 
junction area. Our samples were fabricated on the same 
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chip. Thorough characterization before the actual quan- 
tum measurements revealed that the junctions exhibit a 
very high quality. We showed that we could significantly 
decrease the damping at frequencies relevant for MQT by 
using lumped element inductors, which allowed us to per- 
form our study in the low damping limit. The crossover 
from the thermal to the quantum regime was found to 
have a clear and systematic dependence on junction size, 



which is in perfect agreement with theory. 
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